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Abstract 

We analyze mechanisms and regimes of wave packet spreading in nonlinear 
disordered media. We predict that wave packets can spread in two regimes 
of strong and weak chaos. We discuss resonance probabilities, nonlinear 
diffusion equations, and predict a dynamical crossover from strong to weak 
chaos. The crossover is controlled by the ratio of nonlinear frequency shifts 
and the average eigenvalue spacing of eigenstates of the linear equations 
within one localization volume. We consider generalized models in higher 
lattice dimensions and obtain critical values for the nonlinearity power, the 
dimension, and norm density, which influence possible dynamical outcomes 
in a qualitative way. 
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1. Introduction 

In this paper we will discuss the mechanisms of wave packet spreading 
in nonlinear disordered systems. More specifically, we will consider cases 
when i) the corresponding linear wave equations yield Anderson localization, 
ii) the localization length is bounded from above by a finite value, iii) the 
nonlinearity is compact in real space and therefore does not induce long range 
interactions between eigenstates of the linear equations. 

There are several reasons to analyze such situations. First, wave propa- 
gation in spatially disordered media has been of practical interest since the 
early times of studies of waves. In particular, it became of much practical 
interest for the conductance properties of electrons in semiconductor devices 
more than half a century ago. It was probably these issues which moti- 
vated P. W. Anderson to perform his groundbreaking studies on what is now 
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called Anderson localization [1]. With evolving technology, wave propagation 
became of importance also in photonic and acoustic devices j^yf. Finally, 
recent advances in the control over ultracold atoms in optical potentials made 
it possible to observe Anderson localization there as well |4j. Peter Hanggi 
and collaborators studied properties of wave propagation in disordered me- 
dia by phase space visualization of the underlying dynamical system in high 
dimensional phase spaces (if. 

Second, in many if not all cases wave-wave interactions are of importance, 
or can even be controlled experimentally. Screening effects can reduce the 
long range character of these interactions considerably for electrons. Short 
range interactions also hold for s-wave scattering of atoms. When many 
quantum particles interact, mean field approximations often lead to effective 
nonlinear wave equations. As a result, nonlinear wave equations in disordered 
media become of practical importance. 

Third, there is fundamental interest in understanding, how Anderson lo- 
calization is modified for nonlinear wave equations. All of the above mo- 
tivates the choice of corresponding linear wave equations with finite upper 
bounds on the localization length. Then, the linear equations admit no trans- 
port. Analyzing transport properties of nonlinear disordered wave equations 
allows to observe and characterize the influence of wave-wave interactions on 
Anderson localization in a straightforward way. 

A number of studies was recently devoted to the above subject 0, Q, 0, 
lit, Hoi U, 12, 13, 14 1 . In the present work we will present a detailed analysis 



of the chaotic dynamics which is at the heart of the observed destruction of 
Anderson localization. In particular, we will show that an optional interme- 
diate strong chaos regime of subdiffusive spreading is followed by an even 
slower subdiffusive spreading process in the regime of weak chaos. 



2. Wave equations 

We will use the Hamiltonian of the disordered discrete nonlinear Schrodinger 
equation (DNLS) 

i 

with complex variables xpi, lattice site indices I and nonlinearity strength 
> 0. The random on-site energies e/ are chosen uniformly from the interval 



2 



[ — 2 ' t] ' w ith W denoting the disorder strength. The equations of motion 
are generated by ipi = 8Hd / d(iipf): 

iipi = eitpi + (3\ipi\ 2 ipi - ipi+\ - i/Ji-i . (2) 

Eqs. d2J conserve the energy ([TJ) and the norm S = J2i l^Pil 2 - We note that 
varying the norm of an initial wave packet is strictly equivalent to varying /3. 
Eqs. (JTJ) and (J2]) are derived e. g. when describing two-body interactions in 
ultracold atomic gases on an optical lattice within a mean field approximation 



15|, but also when describing the propagation of light through networks of 
coupled optical waveguides in Kerr media [161 ]. 

Alternatively we also refer to results for the Hamiltonian of the quartic 
Klein-Gordon lattice (KG) 

Uk = £ f + \u\ + \u\ + ±(u w - m)\ (3) 
i 

where ui and pi are respectively the generalized coordinates and momenta, 
and ii are chosen uniformly from the interval [§,§]• The equations of motion 
are ui = —dUx / dui and yield 

ill = -CM - uf + -^( u i+i + u i-i ~ 2m «) • (4) 

Equations @ conserve the energy (J3J). They serve e.g. as simple models 
for the dissipationless dynamics of anharmonic optical lattice vibrations in 



molecular crystals [17|. The energy of an initial state E > serves as a 
control parameter of nonlinearity similar to (3 for the DNLS case. For small 
amplitudes the equations of motion of the KG chain can be approximately 



mapped onto a DNLS model [18] . For the KG model with given parameters 
W and E, the corresponding DNLS model (jTJ) with norm 5 = 1, has a 
nonlinearity parameter /3 w 3WE. The norm density of the DNLS model 
corresponds to the normalized energy density of the KG model. 

The theoretical considerations will be performed within the DNLS frame- 
work. It is straightforward to adapt them to the KG case. 

3. Anderson localization 

For (3 = with ipi = Aiexp(-iXt) Eq. ([T|) is reduced to the linear eigen- 
value problem 

\Ai = eiA t - - A l+1 . (5) 
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The normalized eigenvectors A u i = 1) are the NMs, and the eigen- 

values A„ are the frequencies of the NMs. The width of the eigenfrequency 
spectrum X u of © is A = W + 4 with X v e [-2 - f , 2 + f ] . 

The asymptotic spatial decay of an eigenvector is given by ~ e - l /£( x ») 
where £(A„) is the localization length and ~ 24(4 — \ 2 )/W 2 for weak 

disorder W < 4 [l], [l9j]. The NM participation number p u = ^/J2i^ti * s 
one possible way to quantize the spatial extend (localization volume) of the 
NM. The localization volume V is on average of the order of 3£(0) for weak 
disorder, and tends to V = 1 in the limit of strong disorder. The average 
spacing d of eigenvalues of NMs within the range of a localization volume 
is therefore of the order of d ~ A/V, which becomes d w AH /2 /300 for 
weak disorder. The two scales d < A are expected to determine the packet 
evolution details in the presence of nonlinearity. 

Due to the localized character of the NMs, any localized wave packet with 
size L which is launched into the system for (3 = , will stay localized for 
all times. If L V, then the wave packet will expand into the localization 
volume. This expansion will take a time of the order of r^ n = 2ir/d. If 
instead L > V, no substantial expansion will be observed in real space. We 
remind that Anderson localization is relying on the phase coherence of waves. 
Wave packets which are trapped due to Anderson localization correspond to 
trajectories in phase space evolving on tori, i.e. quasiperiodically in time. 

4. Adding nonlinearity 

The equations of motion of ([2]) in normal mode space read 

(6) 

with the overlap integral 

Iv,V\,V2,Vz ^ ^ A v lA v -^ [ I . (7) 

I 

The variables <f) u determine the complex time-dependent amplitudes of the 
NMs. 

The frequency shift of a single site oscillator induced by the nonlinearity 
is Si = (3\ipi\ 2 - If instead a single mode is excited, its frequency shift can be 
estimated by 5 U = f3\4> v \ 2 /p v . 
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As it follows from (J5]), nonlinearity induces an interaction between NMs. 
Since all NMs are exponentially localized in space, each normal mode is ef- 
fectively coupled to a finite number of neighbouring NMs, i.e. the interaction 
range is finite. However the strength of the coupling is proportional to the 
norm density n = |0| 2 . Let us assume that a wave packet spreads. In the 
course of spreading its norm density will become smaller. Therefore the ef- 
fective coupling strength between NMs decreases as well. At the same time 
the number of excited NMs grows. 

One possible outcome would be: (I) that after some time the coupling 
will be weak enough to be neglected. If neglected, the nonlinear terms are 
removed, the problem is reduced to the linear wave equation, and we obtain 
again Anderson localization. That implies that the trajectory happens to be 
on a quasiperiodic torus. Then it has to be on that torus from the beginning. 
Another possibility is: (II) that spreading continues for all times. That would 
imply that the trajectory evolves not on a quasiperiodic torus, but in some 
chaotic part of phase space. A third possibility is: (III) that the trajectory 
was initially strongly chaotic, but manages in the course of spreading to get 
trapped between denser and denser torus structures in phase space after some 
spreading, leading again to localization as an asymptotic outcome. 

Consider a wave packet with size L and norm density n. Replace it by 
a finite system of size L and norm density n. Such a finite system will be 
in general nonintegrable. Therefore the only possibility to generically obtain 
a quasiperiodic evolution is to be in the regime where the KAM theorem 
holds. Then there is a finite fraction of the available phase space volume 
which is filled with KAM tori. For a given L it is expected that there is a 
critical density tikam{L) below which the KAM regime will hold. We do not 
know this L-dependence. Computational studies may not be very conclusive 
here, since it is hard to distinguish a regime of very weak chaos from a strict 
quasiperiodic one on finite time scales. 

The above first possible outcome (I) (localization) will be realized if the 
packet is launched in a KAM regime. Whether that is possible at all for 
an infinite system is an open issue. The second outcome (II) (spreading) 
implies that we start in a chaotic regime and remain there. Since the packet 
density is reduced and is proportional to its inverse size L at later times, 
this option implies that the critical density tikam{L) decays faster than 
possibly faster than any power of 1/L. The third possibility (III) (asymptotic 
localization) should be observable by some substantial slowing down of the 
spreading process. 
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4-1. The secular normal form 

Let us perform a further transformation <ft u = e~ lX ' A Xu and insert it into 
Eq. ©: 

The right hand side contains oscillating functions with frequencies 

K,u = K + - K 2 - K 3 , V = (z/ 1; I/ 2 , IA5) . (9) 

For certain values of u, V the value becomes exactly zero. These secular 
terms define some slow evolution of OH]). Let us perform an averaging over 
time of all terms in the rhs of (JSJ, leaving therefore only the secular terms. 
The resulting secular normal form equations (SNFE) take the form 

^Xv P ^ ^ Iy,vv\,v\ \Xvi I Xv ■ (10) 

Note that possible missing factors due to index permutations can be absorbed 
into the overlap integrals, and are not of importance for what is following. 
The SNFE can be now solved for any initial condition Xv(t — 0) = f] u and 
yield 

Xu(t)=T) v e-* lvt , n„ = pY.^u^M 2 . (ii) 

Since the norm of every NM is preserved in time for the SNFE, it follows that 
Anderson localization is preserved within the SNFE. The only change one 
obtains is the renormalization of the eigenfrequencies \ v into \ v = \ v + Q u . 
Moreover, the phase coherence of NMs is preserved as well. Any different 
outcome will be therefore due to the nonsecular terms, neglected within the 
SNFE. 

4-2. Measuring properties of wave packets 

We order the NMs in space by increasing value of the center-of-norm co- 
ordinate X v = Yli^ui- We analyze normalized distributions n v > using 
the second moment m 2 = Ylv( v ~ v) 2 n u , which quantifies the wave packet's 
degree of spreading and the participation number P = 1/ Yl v n vi which mea- 
sures the number of the strongest excited sites in n u . Here v = ^2, u yri v . 
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We follow norm density distributions n u = \4>u\ 2 /J2^ I^mI 2 - The second mo- 
ment m2 is sensitive to the distance of the tails of a distribution from the 
center, while the participation number P is a measure of the inhomogeneity 
of the distribution, being insensitive to any spatial correlations. Thus, P 
and 777,2 can be used to quantify the sparseness of a wave packet through the 
compactness index 

p2 

C = — • (12) 
m 2 

A thermalized wave packet yields ( = 3. Distributions with larger gaps 
between equally excited isolated sites attain a compactness index ( < 3. 

4-3. Expected regimes of spreading 

Previous studies suggested the existence of various dynamical regimes 
of spreading of wave packets 0, EU H|- Some of these definitions were 



contradictory. Below we will resolve this. 

Consider a wave packet at t — which has norm density n and size L. If 
f3n > A, a substantial part of the wave packet will be selftrapped 

This 

is due to the above discussed nonlinear frequency shifts, which will tune the 
excited sites immediately out of resonance with the nonexcited neighborhood. 
As a result, discrete breather like structures will be formed, which can persist 
for immensely long times. While selftrapping and discrete breather formation 
are interesting localization phenomena at strong nonlinearity, they are very 
different from Anderson localization since they require the existence of gaps 



in the spectrum of the linear wave equations [20[. If now (3n < A, selftrapping 
is avoided, and the wave packet can start to spread. For L < V and (3 = 0, 
the packet will spread over the localization volume during the time r^ n . 
After that, the new norm density will drop down to n(r^ n ) rj n^. For 
L > V the norm density will not change appreciably up to Tu n , n{ri in ) w 
n. The nonlinear frequency shift (3n(Tu n ) can be now compared with the 
average spacing d. If j3n(Tu n ) > d, all NMs in the packet are resonantly 
interacting with each other. This regime will be coined strong chaos. If 
instead f3n(Ti in ) < d, NMs are weakly interacting with each other. This 
regime will be coined weak chaos. To summarize: 

/3n(Tn n ) < d : weak chaos 
/3n(Tu n ) > d : strong chaos 
(3n > A : selftrapping 
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In terms of the above wave packet characteristics n , L it follows 



f3nL < A : weak chaos 
(3nL > A : strong chaos (13) 
(3n > A : selftrapping 

where L = L for L < V and L = V for L > V. It follows that the regime 
of strong chaos can be observed only if L > 1. For L = 1 we expect only 
two regimes - selftrapping and weak chaos. Furthermore, we obtain that the 
regimes of strong and weak chaos are separated by the quantity (3n = d, i.e. 
the average spacing d is the only characteristic frequency scale here. 

4- 4- Discussion of numerical results 

Let us discuss the above in the light of published computational exper- 



iments. We show results for single site excitations from 11( in FigJT] with 
W = 4, L = 1 and n = 1. For the DNLS model (left plots in FigJU with 
(3 = 4.5 it follows 0n = 4.5. Already at these values selftrapping of a part 
of the wave packet is observed. Therefore P does not grow significantly, 
while the second moment ni2 ~ t a with a k 1/3 (red curves). A part of 
the excitation stays highly localized jof, while another part delocalizes. For 
(3 = 1 selftrapping is avoided since (3n < A. With V ~ 20 and d ~ 0.4 it 
follows that Tun ps 16 and f3n{ri in ) ^ 0.05 d Therefore we expect to ob- 
serve the regime of weak chaos. It is characterized by subdiffusive spreading 
with vri2 ~ t a and P ~ t a l 2 (green curves). For f3 = 0.1 we will remain in 
the regime of weak chaos, however the time scales for observing spreading 
grow. Therefore one finds no visible spreading up to some time which 
increases with further decreasing nonlinearity. For t < both 777,2 and P are 
not changing. However for t > the packet shows visible growth with the 
characteristics of weak chaos (blue curves). The simulation of the equations 
of motion in the absence of nonlinear terms (orange curves) shows Anderson 
localization. Since L = 1 in the above numerical data, strong chaos has not 
been observed. 



Notably, the authors of Ref. [11] also considered single mode excitations 
with total norm £7=1. Using the above terminology, n ps 1/V and L = V 
with W = 4 and therefore again V ~ 20. For the case /3 = 6.5 the authors 
detected a growth of mi which was subdiffusive but faster than t 1//3 . We 
think that these observations are a clear hint towards the realization of strong 
chaos, which should be observable for 5... 10 < < 30... 40 in these cases. 
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Figure 1: (color online) Single site excitations, and P versus time in log-log plots. 
Left plots: DNLS with W = 4, (3 = 0,0.1,1,4.5 [(o), orange; (b), blue; (g) green; (r) red]. 
Right plots: KG with W — 4 and initial energy E = 0.05, 0.4, 1.5 [(b) blue; (g) green; (r) 
red]. The orange curves (o) correspond to the solution of the linear equations of motion, 
where the term uf in Q was absent. The disorder realization is kept unchanged for each 
of the models. Dashed straight lines guide the eye for exponents 1/3 (7712) and 1/6 (P) 
respectively. Insets: the compactness index ( as a function of time in linear-log plots for 
(3 = 1 (DNLS) and E = 0.4 (KG). Adapted from [ll| 
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The time evolution of ( for excitations in the regime of weak chaos is 
shown in the insets of Fig. [TJ As one can see the compactness index oscillates 
around some constant nonzero value both for the DNLS and the KG models. 
This means that the wave packet spreads but does not become more sparse. 
The average value ( of the compactness index over 20 realizations of single 
mode excitations at t = 10 s for the DNLS model with W = 4 and (3 = 5 was 



found to be C = 2.95 ± 0.39 11 



The norm density distribution for the DNLS model at t = 10 8 is plotted 
in Fig. [2j The distribution is characterized by a flat plateau of almost ideally 
thermalized NMs. The width of this plateau is more than an order of mag- 
nitude larger than the localization volume of the linear equations. Therefore 
Anderson localization is destroyed. The plateau is bounded by exponentially 
decreasing tails, with exponents corresponding to the localization length of 
the linear equations. With growing time the plateau widens, drops in height, 
and is pushing the tails to larger distances. Another remarkable feature are 
the huge fluctuations of norm densities in the tails, reaching 4-6 orders of 
magnitude. Such fluctuations are observed even in the case (3 = 0. They 
are due to the fact, that NMs are ordered in space. Neighbouring NMs in 
space may have different eigenfrequencies, and therefore different values of 
their localization length. Tail NMs are excited by the core. The further away 
they are, the weaker the excitation. But within a small tail volume, NMs 
with larger localization length will be more strongly excited than those with 
smaller localization length, hence the large observed fluctuations, which on a 
logarithmic scale are of the order of the relative variation of the localization 
length. The remarkable observation is, that these fluctuations in the tails 
persist for the nonlinear case. Anderson localization is destroyed in the core 
(plateau) of the wave packet due to mode-mode interactions. The tail NMs 
are slaved to the core and excited by it. The interaction between neigh- 
bouring tail NMs is negligible, and the huge fluctuations persist. Therefore, 
Anderson localization is preserved in the tails of the distributions over very 
long times (essentially until the given tail volume becomes a part of the core). 

For single site excitations in the regime of weak chaos the exponent a 
of subdiffusive spreading does not appear to depend on (3 in the case of the 
DNLS model or on the value of E in the case of KG. We find no visible 
dependence of the exponent a on W. Therefore the subdiffusive spreading 
is rather universal and the parameters (3 (or E) and W are only affecting 
the prefactor. Excluding selftrapping, any nonzero nonlinearity appears to 
completely delocalize the wave packet and destroy Anderson localization. 
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Figure 2: Norm density distributions in the NM space at time t = 10 8 for the initial exci- 
tations in the regime of weak chaos of the DNLS model. Left plots: single site excitation 
for W = 4 and /3 = 1. Right plots: single mode excitation for W — 4 and j3 — 5. |<?V| 2 is 
plotted in linear (logarithmic) scale in the upper (lower) plots. The average localization 
volume V fa 20 (shown schematically in the lower plots) is much smaller than the length 
over which the wave packets have spread. Adapted from (Tlj . 
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Fittings were performed by analyzing 20 runs in the regime of weak chaos 
with different disorder realizations. For each realization the exponent a was 
fitted, and then averaged over all computational measurements. We find 



a = 0.33 ± 0.02 for DNLS, and a = 0.33 ± 0.05 for KG [10J, LL3J • Therefore, 
the universal exponent a = 1/3 [lo[ appears to explain the data. 

Another intriguing test was performed on the same disorder realizations 
and single site initial conditions, by additionally dephasing the NMs in a 



random way every hundred time units (11]. In that case, subdiffusion speeds 



up, and m 2 grows as t 1 ^ 2 implying ad ep h = 1/2. This regime of complete 
decoherence of NM phases exactly corresponds to the above anticipated one 
of strong chaos, but here enforced by explicit dephasing. 



5. From strong to weak chaos, from resonances to nonlinear diffu- 
sion 

We can think of two possible mechanisms of wave packet spreading. A 
NM with index /i in a boundary layer of width V in the cold exterior, which 
borders the packet, is either incoherently heated by the packet, or resonantly 
excited by some particular NM from a boundary layer with width V inside 
the packet. 

For heating to work, the packet modes 4> v {t) should contain a part 4>l(t), 
having a continuous frequency spectrum (similar to a white noise), in addition 
to a regular part (fil{t) of pure point frequency spectrum: 

Mt) = + . (14) 

Therefore at least some NMs of the packet should evolve chaotically in time. 
The more the packet spreads, the less the mode amplitudes in the packet 
become. Therefore its dynamics should become more and more regular, im- 
plying lim^oo^ {t)/4>l{t) 0. 

5.1. Strong chaos 

Let us assume that all NMs in the packet are strongly chaotic, and their 
phases can be assumed to be random on the time scales of the observed 
spreading. According to ([6]) the heating of the exterior mode should evolve 
as 20 M « A M M + (3n 3 / 2 f(t) where {f(t)f(f)) = 5(t - f) ensures that f(t) has 
a continuous frequency spectrum. Then the exterior NM increases its norm 
according to |0 M | 2 ~ (3 2 nH. The momentary diffusion rate of the packet is 
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given by the inverse time T it needs to heat the exterior mode up to the 
packet level: D = 1/T ~ j3 2 n 2 . The second moment is of the order of 
m2 ~ l/"^ 2 ; since the packet size is 1/n. The diffusion equation mi ~ Dt 
yields m% ~ f3t 1 / 2 . This agrees very well with the numerical results for 
dephasing in NM space. Moreover, we expect it to hold also without explicit 
dephasing, provided the initial wave packet satisfies the above conditions for 
strong chaos (fT3"j) . First numerical tests show that this is correct 2l|, but 
it contradicts the observations of the numerical data in the regime of weak 
chaos without additional dephasing. Thus, in the regime of weak chaos, not 
all NMs in the packet are chaotic, and dephasing is at best some partial 
outcome. 



5.2. Resonance probability 

Chaos is a combined result of resonances and nonintegrability. Let us 
estimate the number of resonant modes in the packet for the DNLS model. 
Excluding secular interactions, the amplitude of a NM with \(f) u \ 2 = n v is 
modified by a triplet of other modes JL = (fix, fi2, 1^3) m first order in /3 as (jH]) 



(1)1 

V I 



A 



(15) 



where X u ^ = X u + A Ml — A M2 — A M3 . The perturbation approach breaks down, 
and resonances set in, when y/n^ < \<f>v\. Since all considered NMs belong 
to the packet, we assume their norms to be equal to n for what follows. Then 
the resonance condition for a given NM with index v is met if there is at least 
one given triplet of other NMs \x such that 



{3n < R Vjl 



(16) 



If three of the four mode indices are identical, one is left with interacting 
NM pairs. A statistical analysis of the probability of resonant interaction 



was performed in Ref. [10|]. For small values of n (i.e. when the packet 
has spread over many NMs) the main contribution to resonances are due 
to rare multipeak modes jlQ], with peak distances being larger than the 
localization volume. However pair resonances are expected not to contribute 
to the spreading process 22] . When distances between the peaks of multipeak 
modes are larger than the localization volume, the time scale of excitation 
transfer from one peak to another will grow exponentially with the distance. 
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Such processes are too slow in order to be observed in numerical experiments 

3- 

If two or none of the four mode indices are identical, one is left with 
triplets and quadruplets of interacting NMs respectively. In both cases the 
resonance condition f fT6|) can be met at arbitrarily small values of n for NMs 
from one localization volume. 

For a given NM v we define R v ,fi = min^ R v ^. Collecting R u ^ for 
many v and many disorder realizations, we can obtain the probability density 
distribution W(R V ^ ). The probability V for a mode, which is excited to a 
norm n (the average norm density in the packet), to be resonant with at least 
one triplet of other modes at a given value of the interaction parameter (3 is 
therefore given by 

rfin 

V= W(x)dx . (17) 
Jo 



The main result is that W(R u ,jx ->■ 0) ->• C(W) ^ 11]. For the cases 
studied, the constant C drops with increasing disorder strength W . This 
result of nonzero values of C is not contradicting the fact of level repulsion 
of neighbouring NMs, since triplet and quadruplet combinations of NM fre- 
quencies can yield practically zero values of with finite distances between 
the eigenfrequencies. 

For the case of strong disorder (W ^> 1) the localization volume tends 
to one, and quadruplet resonances are rare. Excluding also pair resonances 
for the above reasons, we are left with triplet resonances. A given mode 
may yield a triplet resonance with its two nearest neighbours to the left and 
right. Replacing the overlap integrals by some characteristic average, and 
assuming that the three participating modes have essentially uncorrelated 
eigenfrequencies, it follows that 

W(R)nC (l-9^L\ . (is) 



Due to the nonnegativity of V it would follow that V = for R > 3/C. In 
reality we expect an exponential tail for large R. As a simple approximation, 
we may instead use 

W(R) w Ce~ CR (19) 

which in turn can be expected to hold also for the case of weak disorder. It 
leads to the approximative result 

V = l- e~ ci3n . (20) 
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For /3n — > it follows 

V « . (21) 
Therefore the probability for a mode in the packet to be resonant is propor- 



tional to Cf3n in the limit of small n [10|,lll|- However, on average the number 



of resonant modes in the packet is proportional to the product of V and the 
total number of modes in the packet. Since the total number is proportional 
to 1/n, the the average number of resonant modes in a packet is constant, 
proportional to C/3, and their fraction within the packet is ~ c(3n (num. 
Since packet mode amplitudes fluctuate in general, averaging is meant both 
over the packet, and over suitably long time windows. A detailed numerical 
analysis of the statistical properties of resonances and related issues is in 



preparation [23 



Finally we consider the process of resonant excitation of an exterior mode 
by a mode from the packet. The number of packet modes in a layer of the 
width of the localization volume at the edge, which are resonant with a 
cold exterior mode, will be proportional to /3n. After long enough spreading 
j3n <C 1. On average there will be no mode inside the packet, which could 
efficiently resonate with an exterior mode. Resonant growth can be excluded 



10|, lU|. Thus, a wave packet is trapped at its edges, and stays localized until 
the interior of the wave packet decoheres (thermalizes) . On these (growing) 
time scales, the packet will be finally able to incoherently excite the exterior 
and to extend its size. 



5.3. A conjecture leading to the correct asymptotics 

We assume, that the continuous frequency part of the dynamics of a 
packet mode is V((3n). It follows that 4>l(t) / '4>l(t) ~ V((3n). As expected 
initially, the chaotic part in the dynamics of packet modes becomes weaker 
the more the packet spreads, and the packet dynamics becomes more and 
more regular in the limit of large times. Therefore the chaotic component 
is conjectured to be a small parameter 4> c u {t) <C Expanding the term 

|0v| 2 0i/ to first order in 4>l(t), the heating of the exterior mode should evolve 
according to « A M M + /3n 3 / 2 V((3n)f(t). It follows |0 M | 2 ~ (3 2 n 3 (V((3n))H, 
and the rate 

D = 1/T ~ /3 2 n 2 (V(/3n)) 2 . (22) 
With (!T9l).(l20l) and ~ 1/n 2 the diffusion equation mi ~ Dt yields 

\ ~ 0(1 - e- c/3n )t 1/2 . (23) 
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The solution of this equation yields a crossover from subdiffusive spreading 
in the regime of strong chaos to subdiffusive spreading in the regime of weak 
chaos: 

1^2 ~ (/3 2 t) 1 ^ 2 , strong chaos , Cf3n > 1 , 
m 2 ~ C 2/3 /3 4/3 t 1/3 , weak chaos , C/3n < 1 , 

5.4- The crossover from strong chaos to weak chaos 

According to (ITT)) the probability of resonance for a packet NM will be 
practically equal to one, if (3n is sufficiently larger than 1/C. Such a situation 
can be generated for packets with large enough fin, and should yield spread- 
ing, provided one avoids selftrapping f3n < A + W P, [llj . This spreading will 
be different from the asymptotic behaviour discussed above over potentially 
large time scales. 

Let us use as an example W = 4 and (3 = 1, with the constant C ~ 6.2 



Single site excitations with norm S = 1 lead after very short times 
to a spreading of the excitation into the localization volume of the linear 
wave equations, which is of the order of 10-20. The attained norm density is 
therefore of the order of n < 0.1. The observed spreading is the asymptotic 
one since V ~ Cf3n. However, if we choose a packet size L to be of the order 
of the localization volume, and the norm density n of the order of n = 1, 
initially V ~ 1. Thus every mode in the packet will be resonant, and the 
condition for strong chaos should hold. At the same time (3n = 1 is far 
below the selftrapping threshold 4 + W = 8. For strong chaos we derived 
?77, 2 t 1 / 2 . With spreading continuing, the norm density in the packet will 
decrease, and eventually (3n < 1/C. Then there will be a crossover from 
strong chaos to weak chaos, and m 2 ~ t 1 ^ 3 for larger times. This crossover 
happens on logarithmic time scales, and it will be not easy to confirm it 



numerically [21 



In Figj3] we show the resulting time dependence of rri2 on t from (|23|) in 
a log-log plot, where we used (3 = 1, C = 6.2, L = 20 and n(t = 10 2 ) = 1. 
With x = log 10 (t) and y = /o^io(^ 2 ) it is straightforward to calculate the 
zero of the third derivative d 3 y/dx 3 = to obtain the crossover position: 

C(3n c w 1.86 . (24) 

Therefore the only characteristic frequency scale here is 1/C. From the above 
discussion of the different spreading regimes (IT3)) it follows, that this scale is 
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Figure 3: m 2 (t) in a log-log plot according to (|23|) (black solid line). Dashed lines - power 
laws for strong and weak chaos. 
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corresponding to the average spacing d: 



i « d . (25) 



Then 



C~^ 2 ,W<4, (26) 

C » , W » 4 . (27) 

Our results can be used to predict the critical value of the norm density n c 
at which the crossover should take place. For W = 4 and = 1 it follows 
n c « 0.3. 

5. 5. Scaling of nonlinear diffusion equations 

With the above results on the diffusion coefficient, we may consider a 
set of nonlinear diffusion equations for the norm density distributions in NM 
space. For simplicity we replace the discrete NM indices by a continuous 
variable: 

dn d „ . .dn . . 

at = • (28) 

In particular, we are interested in cases D(n) ~ n K . With a single scaling 
assumption n(u,t/a) = bn(cv, t), and using the conservation of the total 
norm and n{y — > ±oo,t) — > we obtain b = c = a l '^ K+2 '. Then the second 
moment m 2 will grow in time according to 

Mt) = (f) m 2 (t ) , a = . (29) 

Notably an explicit self-similar solution was calculated by Tuck in 1976 
24j which has the following spatial profile: 

2 \ V« 

2(k + 2) 



^)=(5-^^t) • (30) 



Here i? is an integration constant (see also (25j). 

For the case of strong chaos k = 2 and a = 1/2 in accord with the above 
results. Also for k = 4 we obtain a = 1/3 again in agreement with the above 
results. 



18 



6. Generalizations 



Let us consider D- dimensional lattices with nonlinearity order a > 0: 

Here I denotes an D- dimensional lattice vector with integer components, and 
m G D(l) defines its set of nearest neighbour lattice sites. We assume that 
(a) all NMs are spatially localized (which can be obtained for strong enough 
disorder W), (b) the property W(x — > 0) — > const ^ holds, and (c) the 
probability of resonances on the edge surface of a wave packet is tending to 
zero during the spreading process. A wavepacket with average norm n per 
excited mode has a second moment m 2 ~ \jn 2 l r> . The nonlinear frequency 
shift is proportional to /3?W 2 . The typical localization volume of a NM is 
still denoted by V, and the average spacing by d. 

Consider a wave packet with norm density n and volume L < V . A 
straightforward generalization of the expected regimes of spreading leads to 
the following: 

cr/2 



/ L\ a/ 

f3n a/2 ( — J V < A : weak chaos 
/ L \ a/2 

j3n u l 2 [y\ V > A : strong chaos 

f3n a / 2 > A : selftrapping 

The regime of strong chaos, which is located between selftrapping and weak 
chaos, can be observed only if 

L>L C = V 1 - 2 ^, n>n c = ^[^j l(T . (32) 

For a = 2 we need L > 1, for a — > oo we need L > V, and for a < 2 we 
need L > 1. Thus the regime of strong chaos can be observed e.g. in a 
one-dimensional system with a single site excitation and a < 2. 

If the wave packet size L > V then the conditions for observing different 
regimes simplify to 

f3n°~/ 2 < d : weak chaos 
f3n°~/ 2 > d : strong chaos 
/3n a/2 > A : selftrapping 
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The regime of strong chaos can be observed if 



d 



2/<7 



n>nc= \j) ■ (33) 

Similar to the above we obtain a diffusion coefficient 

D ~ f3 2 n a (V(/3n a/2 )) 2 . (34) 
In both regimes of strong and weak chaos the spreading is subdiffusive (lo| : 

m 2 ~ (/3 t) 2+<T - D , strong chaos , (35) 
m 2 ~ (/3 4 ^) 1+CT - D , weak chaos . (36) 

Let us calculate the number of resonances in the wave packet volume 
(Nrv) and on its surface (Nrs) in the regime of weak chaos: 

N RV ~ Pn°' 2 - 1 , N RS ~ Pn* 2 ^ 1 . (37) 

We find that there is a critical value of nonlinearity power a c = 2 such that 
the number of volume resonances grows for a < o c with time, drops for 
a > a c and stays constant for a = a c . Therefore subdiffusive spreading will 
be more effective for o < a c . 

We also find that the number of surface resonances will grow with time 

for 

D> D c = 1 —— , a < 2 . (38) 

Therefore, for these cases, the wave packet surface will not stay compact. 
Instead surface resonances will lead to a resonant leakage of excitations into 
the exterior. This process will increase the surface area, and therefore lead 
to even more surface resonances, which again increase the surface area, and 
so on. The wave packet will fragmentize, perhaps get a fractal-like structure, 
and lower its compactness index. The spreading of the wave packet will speed 
up, but will not anymore be due to pure incoherent transfer, instead it will 
become a complicated mixture of incoherent and coherent transfer processes. 

Mulansky computed spreading exponents for single site excitations with 
/3 = 1, W = 4, L = 1, D = 1 n = 1 and a = 1,2,4,6 [13j. Since for 
a = 2, 4, 6 strong chaos is avoided, the fitting of the dependence m s (t) with a 
single power law is reasonable. The corresponding fitted exponents 0.31±0.04 
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(a = 2), 0.18 ± 0.04 (a = 4) and 0.14 ± 0.05 (tr = 6) agree well with the 
predicted weak chaos result 1/3,1/5,1/7 from (1551) . For a = 1 the initial 
condition is launched in the regime of strong chaos. A single power law fit 
will therefore not be reasonable. Since the outcome is a mixture of first strong 
and later possibly weak chaos, the fitted exponent should be a number which 
is located between the two theoretical values 1/2 and 2/3. Indeed, Mulansky 
reports a number 0.56 ± 0.04 confirming our prediction. Veksler et al [12| 
considered short time evolutions of single site excitations (up to t = 10 3 ). 
While the time window may happen to be too short for conclusive results, 
the observed increase of fitted exponents with increasing (3 for o < 2 is 
possibly also influenced by the crossover from weak to strong chaos. Note 
that Skokos et al 26j performed a more detailed analysis for the KG lattice, 
which confirm many the above results. 



7. Discussion and conclusions 

Let us summarize the findings. If the strength of nonlinearity is large 
enough, a wave packet (or at least an appreciable part of it) is selftrapped 
due to the finite bounds for the spectrum of the linear equation. If the 
nonlinearity is weak enough so as to avoid selftrapping, two possible outcomes 
are predicted, which now depend also on the volume L of the packet. If 
L > L c and n > n c , the NMs in the packet will be all resonant, strongly 
interacting with each other and quickly dephase. That leads to a regime of 
strong chaos. As time grows, the norm density n drops below n c , and the 
spreading continues in the regime of weak chaos. If however either L < L c or 
n < n c , strong chaos is avoided, and the packet will spread in the regime of 
weak chaos. Lowering or n further will keep the spreading in the regime of 
weak chaos, but time scales of subdiffusion will grow, and the process will not 
be observable on the finite time window currently accessible by computational 
experiments. The above holds if D < D c which implies that Anderson 
localization is preserved in the tails and destroyed in the wave packet core. 
In other words, the time scales for destroying Anderson localization in the 
tails are much larger than the time scales which lead to a thermalization of 
the core and the corresponding subdiffusive spreading. In order to observe 
the crossover from strong to weak chaos, one has to carefully choose the 
system parameters. In particular, it is desirable to make the crossover region 
more narrow. 

If D > D c then the spreading process will be different from the above 
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predictions, because resonant interaction in the surface and the tails of the 
wave packet will destroy Anderson localization as well. The spreading will 
presumably stay subdiffusive. But we do not know currently how to estimate 
and characterize the details of this process. 

Our results rely on a conjecture of the dependence of a diffusion coeffi- 
cient on the probability of resonances. Future investigations may consider 
the connection between this conjecture and the dependence of Lyapunov co- 
efficients, relaxation times of correlation functions, and detrapping times on 
the system parameters. 
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